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SECTION I 


INTRODUCTION 


Advanced fiber reinforced composites such as boron-epoxy, graphite- 
epoxy, etc. are being used increasingly in flight vehicle structures. 
In order to make the most efficient use of these composites, it is 
necessary to have a reliable design procedure. In References 1 and 2, 
an optimization method was presented for the minimum weight design of 
composite structures with stress and displacement constraints. The 
concepts developed in these two references for the stress constraint 
design form the basis of a computer program to Optimize Composite (OPTCOMP) 
structures discussed in this report. 
j The OPTCOMP program can be used to design a minimum weight structure 
or to analyze a given structure which is idealized with finite elements. 
This is an experimental program and as such it includes different recur- 
rence relations and strength criteria for experimentation. A metal 
structure or a combination of metal and composite structure can be 
designed by oniper definition of the elastic properties of the elements. 
The basic relations of the finite element method, the optimality 
criteria, the recurrence relations and the various strength criteria are 


discussed in Section II. The buckling equations and the transformations 


for the elastic constant matrix are also given in Section II. Section III 
contains a general description of the computer program. The input and 
output instructions are included in Sections IV and V. Section VI dis- 
cusses the sample problem. The Fortran listing of the program and the 


' output of the sample problem are presented in an Appendix. 
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SECTION II 


ANALYTICAL FORMULATION 


The optimization of a structure is an iterative procedure. The two 
main steps involved in the algorithm are 1) to analyze the structure to 
determine the response of the structure to the applied loads and 2) to 
modify the design variables by using a recurrence relation. These two 
steps are repeated for a specific number of iterations. The recurrence 
relations used to modify the design variables are based on the selected 
optimality criteria. The optimality criteria used to derive the recurrence 
relation may be theoretical or intuitive. In the OPTCOMP program, the 
response of the structure is obtained by a finite element analysis of the 
structure. The geometry of the structure is fixed. Hence, the design vari- 
ables are the thicknesses of the homogeneous plate elements, the thicknesses 
of the layers in the composite elements, and the cross-sectional areas 
of the bars. As used throughout this report, the word "layer" identifies 
the total quantity of all the plies or laminae which are oriented in a 
specified direction within a composite laminate. Since out-of-plane 
bending effects are not considered in this program, the concept of layers 
may be applied. However, in actual laminates, these plies will be dis- 
persed throughout the thickness. In the composite elements, the fiber 
directions to be used are specified. 

Le Equations of Analysis 

Since the finite element method of analysis is used in OPTCOMP, 

the basic equations are given here for ready reference. In this method 


the structure is discretized into a group of elements. Let the number 


of elements in the structure be m and the number of layers in the compo- 
site elements be n. The force displacement relation for the whole 


structure is given by 
[K] {r} = {R} (1) 


where [K] is the stiffness matrix of the structure, {r} is the displace- 
ment vector and {R} is the force vector. The force displacement relation 


for the jth layer of the ith element is given by 


[k],, {r}, = {fh (2) 


where [ks is the stiffness matrix of the jth layer of the ith element, 
(r}, are the ith element nodal displacements and {f},, are the nodal 
forces of the jth layer of the ith element. The strain energy ei, of 


J 
the jth layer of the ith element is given by 


cd 
ei = {r}5 (kl, {r}5 (3) 


The element strains fe}, in the ith element are given by 

te}, = [c]} (r}, (4) 
where [c] is the strain-displacement matrix. The stresses fo}, in the 
jth layer of the ith element are given by 

fo}, = (Al, fe}, (5) 


where [A], 5 is the matrix of elastic constants of the jth layer of the 
ith element. 
y a Strength Criteria 

The state of stress in a layer of an element has to satisfy a 


certain specified strength criteria in order to avoid failure. There 


age 5 yao 


are several strength criteria proposed by different authors (For a detailed 
discussion, see Reference 3.),. Most of these are explicit or implicit functions 
of the stresses or strains and the corresponding maximum allowable stresses 

or strains. 

For an orthotropic material, there are basically five allowables 
required for each layer of a plane stress problem. These are the tension 
and compression maximum allowables in the fiber direction and in the trans- 
verse direction, and the maximum shear allowable. 

In the OPTCOMP computer program four strength criteria are included. 
The user of the program can select the one he prefers. The four criteria 
for the generalized plane stress condition are given below. 

a. Maximum Stress Criteria 

This criteria states that the failure of the material would occur 


if any one of the following conditions are not satisfied. 


vo ~ os] 
O, S % (6) 
“It, "33 


where o and 55 are the stresses in the fiber direction and the one normal 


to it and 19 is the shear stress. % and J5 are the maximum allowable 


stresses in the fiber and transverse directions, and 19 is the maximum 
allowable shear stress. 
b. Maximum Strain Criteria 


In this criteria strains instead of stresses are checked. The 


criteria can be expressed as 


; = ] 
Eo S Ey (7) 
P49 SS y9 


where ey and €, are the strains in the fiber direction and the one normal 
to it and 19 is the shear strain. aT and E, are the maximum allowable 


strains in the fiber and transverse directions, and me is the maximum 
allowable shear strain. The layer is considered to have failed if any one 
of the conditions specified in Eq. 7 is not satisfied. 

In the maximum stress and strain criteria the individual failure 
modes are independent. 

Cs Hill's Criteria Modified by Tsai 


b According to this criteria the material is assumed not to have 


failed if 


x 

eae \ 2 

[ o, 2 | #4 \ ie o> 19 < 1 (8) 
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| a2 
oy ay *t / 12 
The definitions of op Tos etc. are the same as in the maximum stress 
criteria above. 
d. Norris Criteria 


This criteria can be written as 


This criteria is similar to Eq. 8 except for the third term. 


It should be pointed out that different versions of the above four 
criteria can be used by a proper choice of the allowables. For example, 
if the transverse stress is not to be considered critical as in the case 
of a fiber controlled laminate, 0, can be assigned a high value so that 


3, will not become critical. Any special type of criteria can be easily 


introduced into the program by writing a new subroutine. 
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3. Recurrence Relations 

In this section the recurrence relations incorporated in the OPTCOMP 
program are given. For the stress constraint problem, when the maximum 
allowable stress is not the same in all the elements of a structure, 
there is no rigorous method based on optimality criteria for designing 
a minimum weight structure. This is particularly true when the structure 
has plate elements. If the structure consists of bar elements only, the 
problem can be correctly solved by replacing the stress constraints by 
equivalent displacement constraints and treating it as a multiple dis- 
placement constraint problem (Reference 2). However, this procedure of 
treating stress constraints by equivalent displacement constraints is 
not feasible in a structure where there are elements connecting more 
than two nodes. The recurrence relations proposed here for resizing the 
elements are based on optimality criteria which are not rigorous in the 
mathematical sense, but they are found to give near optimum weight 
designs for large structures in a very efficient way. The recurrence 
relations included in OPTCOMP are either based on the strain energy in 
an element or the state of stress in an element. Those depending on 
the state of stress correspond to the fully stressed iterative procedure. 
The user of the program may select any one of them depending on his 
preference. 

The optimality criteria for the generalized stiffness requirement 
can be stated as, "the optimum structure is the one in which the ratio 
of the average strain energy density to the mass density is the same in 


all elements" (Reference 2). This criteria can be written as 


ae 
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e 
ee a (10) 
Pay 
i=l, --m 
j= ly na 
where Ais the Lagrangian parameter, P45 is the mass density of the jth 
layer of the ith element, and es is the strain energy density given by 
— = 
eu 7 gt (11) 
ij 
Here Vay is the volume of the element defined as 
¥ 44 = tay 1; (12) 
where Cs is the thickness of the jth layer of the ith element and 1 is 
the surface area of the ith element. The design variable ts can be 
written as 
tay = \ ee (13) 
where @ is the relative thickness (normalized to the maximum thickness) 


ij 


of the jth layer of the ith element and A is the scaling parameter. For 
homogeneous elements in the structure j = 1. Introducing the scaling 


parameter A in Eq. 1 gives 


hot] teh ee AR (14) 
or 

[K') {r'} = {R} (15) 
where 

(r} = 2 (r'} (16) 


In Eq. 15 [K'] is the stiffness matrix for the whole structure 


obtained by using the relative design vector a Introducing the scaling 


» ie 


parameter into Eqs. 2 through 5, the relations between the actual quan- 


tities and the relative quantities at element level can be expressed 


SE OR, titra i 7 . eer aeinetenmenatneeesie ee 


—- 


——— 
Le —— 
picecisinlhs SACI AARC AO ES a On — 
t 
‘ 


as 
= ' C17 
(k],, A(k Vi, ) 
{x}, = + {r'} (18) 
uf A i 
ae a re 19 
fe}, rn {e i (19) 
{o},, = + {0} (20) 
ij A ij 
where the primed quantities are the relative values. 
Introducing Eqs. 11, 17 and 18 in Eq. 10 gives 
e'. 
1-4 — (21) 
oes 
' t ' ' 
where ate {r yi [k'],, {r hy (22) 
ij ay, 1 


Multiplying both sides of Eq. 21 by a and taking the square root 


gives 


Zz 
eae = Ba — 


where B is a constant. Eq. 23 can be rewritten in an iterative form as 


at ly 
=. (24) 
ij Vv 


D 


ay j B44) 


Pe 
| 
Vv 
where v+l and v refer to the cycles of iteration. In Eq. 24, Bis a 
constant and generally it is not essential to evaluate it, since only 
the normalized relative design vector ve is used to analyze the struc- 


ture. In the case of multiple loading conditions, a sum of relative strain 


energy densities due to all the loading conditions is used in Eq. 24. 


we. 
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In order to take into consideration the effect of different allowable 


stresses, the optimality criteria defined by Eq. 10 may be modified as 


L= } “2. (25) 
Diy 


where 


a Fa ee | (26) 
(0) Evsy 


where o is the maximum allowable stress and ey is the maximum allowable 
strain in the fiber direction of the jth layer of the ith element. 
Following the procedure similar to the one used to derive Eq. 24, the 


recurrence relation corresponding to the optimality criteria defined 


by Eq. 25 can be written as 


O,, = Bla | ey , (27) 
( tj v+l1 ( 15), Cy 
Vv 
where 
iT ie ae sed C203 
(Oy &y)45 


In Eq. 28 the prime quantities are the relative values. In the case of 
multiple loading conditions et is calculated for each loading condition 


and the sum is used in Eq. 27. 
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The fully stressed iterative procedure is based on the intuition 

that a structure will have minimum weight if the stress in each member 

is equal to the maximum allowable stress for that member, under at 

least one loading condition. In this procedure the resizing of an 

element is done by multiplying the design variable by the ratio of the 
maximum actual stress to the corresponding maximum allowable stress for 

that member. In the case of a plate element since there will be three 
stresses (two normal and one shear), the effective stress ratio can be used. 
Hence, in the case of Hill's criteria modified by Tsai (Eq. 8) the 


recurrence relation can be written as 


“1 
‘oy = (a,. se iad (29) 
rien ( 15), loos : 
where 
2 2 2 
ce o! Oo. oO: 0 
‘ Pe | ioe | 2 Pe hee SOE 8 _12 
© ogg (Ht) = = ay (- ~2 oe : (30) 
"y 2 1 12 | | 
max 


Similarly, to fully stress the design with Norris criteria (Eq. 9), the 


recurrence relation is 


(a)? (1), feel” *) 


where 


(N) 


' 
eff 


} 1\2 fe © 2 

(<i 7 Rio 1 ra o3| + {| °12 | (32) 
\ 
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nax 


1l 


er 
neath nese tia A! ——— ce ta Ney Aa eC Tt tii 
‘ 


In using the recurrence relation to resize a bar element or a shear 
panel, only the appropriate component of the stress vector is used. 
4. Buckling Equations 

The equations used to modify the layer thickness of a plate element 
to prevent local buckling are given in this section (Reference 4). For 
a simply supported orthotropic plate where the length to width ratio is 


large (>>1), 


2 
pe rae ee (33) 
xcr x Z 
w 
and 
x Diy 
N 2c (34) 
} xycr xy we 
i 
; where N and N are the critical compressive load (parallel to 0° 
xer xycer 


fiber) and critical shear load respectively. In Eqs. 33 and 34, 


K.. = 2/8 +a (35) 
and 
a (2.29 + a ph * cor Y (36) 
' xy 8 — 
4 ~9 ap 1/2 of 
Rey = (4.742 + 0.216 6 + 0.380 6°) (28) rer oO < 2. (37) 
where 
2(D © Zp.) 
iL? 
az ate PB. (38) 
ll 
D 
22 
s=— (39) 
et! 
ss og’? 
Ss 
a (40) 
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In Eqs. 33 and 34, w is the effective width of the plate and Dip 
Dio» etc. are the elements of the bending stiffness matrix. When the 
plate is subjected to combined compressive and shear load, the stability 


criterion is 


2 
"x + Oxy <2 (41) 
Nxer xycr 


where ny and Nyy are the actual compressive and shear force per unit 
length applied to the laminate. If the elastic properties are smeared 


One can assume that 


mee: Sas (42) 


where t is the total thickness of the plate. Substituting Eq. 42 in 


Eqs. 33 and 34 gives 


and 


where 


and 


ee 


ee — 


Substituting Eqs. 43 and 44 in Eq. 41 gives 


n 
re a ee | (47) 


: ——— 3 sa 
Ky t Kyy t 


Hence, knowing ny» ty? Ky and Fi Eq. 47 can be solved for t to 
satisfy the stability criteria in the equality sense. In the OPTCOMP 
program during each iteration in the second stage, the required 't"' is 
evaluated and the thickness of each layer is increased proportional to 
the specified fractions. The influence of boundary conditions can be 
taken into consideration by specifying the proper value of the width, w, 
of the plate. 
oP Transformations 

In this section the transformations used to define the fiber 
direction and the elastic constants in the local coordinate system are 
given. The fiber direction for 0° fibers can be defined in OPTCOMP either 
by giving the angle 6 as shown in Figure 1(a) separately for each element 
or by specifying the projection of the 0° fibers on a defined plane (Figure 
1(b)). If the latter is used, the projected direction will be the same for 
all 0° fibers in all the elements. In Figure 1(a) (e,> eo» e) is the 
local coordinate system used to define the position of the quadrilateral 
ABCD. Let AF be the unit vector parallel to the 0° fibers in the element. 
On the plane PQRS in Figure 1(b), A,F, is the projection of AF. The plane 


ae 


PQRS is defined by ALF) and A,G) which are perpendicular to each other. 


n is the normal vector of unit length to the plane PQRS. 
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Consider a plane passing through the vectors AVF, n and AF. (See 


Figure l(c)). Then 


At, + FOF = Ah = AF (48) 


and 


= 
" 
= 

>| 


(49) 


where H is a scalar quantity. The relation given in Eq. 49 is correct 


since the unit vector n is parallel to FYF. Taking the scalar product 


of Eq. 48 with e. and using Eq. 49, one can write 


H = - A. - 1 (50) 
e n 
n 
Thus, knowing the direction of the projected vector ALF, Eqs. 48 


through 50 can be used to define the vector AF. In the input data of 


the OPTCOMP program the vectors AFL and A\G) are defined by giving the 


coordinates of points Ay> Fy and G)- 


In Figure 2, (e, a, ) is the local coordinate system of the 


element. The unit vector parallel to the fiber direction is ee e, is 


perpendicular to e, such that e. =e 


x e . The two coordinate 
f 8 


f 


systems are related by 


where fa bas etc. are the direction cosines. Let {co} and {ce} be the 
stresses and strains in the local coordinate system, and {co} and {€} be 
the same quantities in the fiber coordinate system (ee, ae €.): The 


& 
stresses and the strains in the two coordinate systems are related by 


{o} = [L,] {o} (52) 
{ce} = (L,] fe} (53) 
where r + 
»2 » 2 
i fi 
Mage 12 2ba1 bap 
oe 2 2 ) 
i] | oo bo9 2koo boy 
£ D £ L V4 e e 
| Sr oe ee ee ee ey * ay “ee 
and 
[ 92 2 
| 1l 12 3 Ca Wy: 
| 
oe 2 
(L4] | 21 22 22 21 
cas et ie ee 12 “21 + “na 


-1 t 
L = L 
[15] [5] 
and 
-1 t 
[L5] = [1] 
The stress-strain relations are given by 
{fo} = fA} {6} 
16 


and 


{o} = [A] fe} (59) 
where [A] and [A] are the elastic constant matrices in the two coordinate 


systems. Using Eqs. 52 and 53, Eq. 59 can be written as 
[1,1] fo} = [al [Ly] fe) (60) 
Multiplying both sides of Eq. 60 by 1 ie gives 


ist 8 3 TR [1] {e) (61) 


1 


or 
{o} = (1,1 {Al [L] {e} (62) 
Comparing Eas. 58 and 62 we can write 


A] = [L,]° [A] (L 
(A) = 4,1" (8) (4,1 (63) 
Eq. 63 can be used to transform the elastic constant matrix from one 


coordinate system to the other. 
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SECTION III 


DESCRIPTION OF THE PROGRAM 


The flow diagram showing the basic steps of the optimization proce- 
dure used in OPTCOMP is given in Figure ?. In the first stage, the 
structure is optimized to satisfy the stress constraints alone. The 
design at the end of this stage will have minimum weight, and it will 
satisfy the strength criteria in all the elements. In the second stage, 
the plate elements which are critical in buckling are increased in thick- 
ness to prevent local buckling. At the end of the second stage, all 
elements will satisfy the strength criteria and no plate elements will 
fail due to local buckling. The initial thicknesses of the elements for 
optimization can be included in the input data or the program will assign 
0.1l-inch thickness for plate elements with equal percentage of laminas in 
all fiber directions, and 1.0 square inch cross-sectional area for bar 
elements. 

The four types of elements included in the program are 1) a constant 
strain triangle, 2) a quadrilateral constructed from four constant strain 
triangles, 3) a shear panel, and 4) a bar. The sequence of node numbers 
used to define the elements is given in Figure 4. The triangle and the 
quadrilateral may be layered composite elements. In defining the material 
axis, either the direction of the 0° fiber orientation is specified by 
angle 6 (Figure 1(a)) defined by the variable ZANGLE in the input Card 
Set 11 for each element, or the projection of the 0° fibers is given on 
a certain plane (See Figure 1(b)). This projection will be the same for 
all the elements. The projection and the plane is defined by the coor- 


dinates of the three points in the input Card Set 9. The fiber 
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directions in other layers are given in the input Card Set 10, with 


respect to the 0° direction. The material properties are given for all 
layers in an element as a set, which is designated by one material 
number. For example, an element with all layers of graphite epoxy may 
be assigned material number 1, an element with one layer of boron-epoxy 
and the remaining layers of graphite epoxy may be assigned material 
number 2, and an element with two layers of graphite-epoxy and the 
remaining layers of glass-epoxy may be assigned material number 3. If 
there are any metal elements with isotropic properties they can be 
assigned yet a different material number. The total number of materials 
to be used for quadrilaterals and triangles is defined by the variable 

I NMAT2 in the input Card Set 3, and the elastic properties, thickness of 

' lamina, minimum number of lamina in each layer are read in the input 
Card Set 7. In the case of shear panels and bar elements the number of 
materials used is defined by NMAT1 in Card Set 3 and the elastic proper- 
ties and the maximum allowable stresses are read in Card Set 6. It should 
be pointed out here that the material numbers assigned for the quadri-~ 


laterals and the triangles may be the same as the one used for the shear 


panels and bars, but their meanings are different for the two cases. 


The program has provision for four widely used strength criteria 


for composite structures. These are: 1) maximum allowable stress, 2) | 


maximum allowable strain, 3) Hill's criteria modified by Tsai, and 4) 


Norris criteria. These criteria are given in Section II(2). The user 
can select any one of these criteria depending on his preference by 


specifying the proper number for NCRTIA in the input Card Set 3, and 
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for NEF in the input Card Set 4. All these criteria are functions of 
the maximum normal allowable stress or strain in the fiber direction 

and in the transverse direction, and the maximum shear allowable. ‘The 
maximum normal allowable stresses may have different magnitudes depen- 
ding on whether they are tensile or compressive stresses. As such the 
total number of maximum allowable stresses specified is five. In case 
the user wants to relax the restriction on any of these allowables, he 
can specify a very high value. Selection of these allowables is very 
important since the weight of the structure is directly proportional to 
the scaling parameter which is determined on the basis of satisfying the 
strength criteria for each element. This data is included in the input 
Card Set 8 for each layer of the material set. 

There are four recurrence relations included in the program (See 
Section II(3)). Two relations are based on the strain energy density 
in an element and the remaining are based on the state of stress. The 
user may select any one of these by assigning the proper value to NENG 
in the input Card Set 3. If the user is not familiar with the implications 
of using the different recurrence relations, it is recommended that he use 
NENG = 1. This will produce good results for a general structure. 

In the input Card Set 3, LSTCCL designates the total number of 
iterations for the stress constraint design. Out of this total, for 
LLSTCL iterations, the design variables are modified by using the recur- 
rence relation. In the remaining iterations, the thickness of each element 
is modified to make the scaling parameter A (See Eqs. 13, 19, 20) nearly 
the same for all the elements. During this second phase of optimization, 


the percentage of laminas in the different fiber directions in the composite 
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elements remains unchanged. This process of jwalizing the scaling 
parameter sometimes reduces the weight of the structure depending upon 
the type of structure, the variation in the magnitude of the different 
maximum allowables, and the selected strength criteria. At the end of 
LSTCCL iterations, the details of the minimum weight design are printed 
out. In the second stage for LNSB iterations, the thicknesses of the 
quadrilateral and triangular elements are modified to prevent local 
buckling of the panels (See Section II1(4)). These equations are based 

on the approximation of smeared orthotropic properties. The effect of 
boundary conditions can be taken into consideration by specifying the 
proper value of the effective width of the plate (AWIDE(1), BWIDE(I) in 
the input Card Set 11). It should be noted that the effective width is 
not the width of the element, but the distance between the supports 
multiplied by the proper fraction to account for boundary conditions 
other than simply supported. The user has the option to indicate which 
layer thicknesses to increase and in what proportion to satisfy the 
stability criteria. This is achieved by giving proper values to PECT 

in the input Card Set 5. In the case of a 0°, 90°, + 45° laminate, it 
has been found that it is more beneficial to increase the + 45° layers 
than to increase the thickness of the whole laminate. The initial design 
for the second stage of optimization is the minimum weight design obtained 
in the first stage. In the second stage, since the thicknesses of the 
elements are increased, the weight of the structure goes up with each 
additional iteration. When all the plate elemen':s satisfy the stability 
criteria, the weight of the structure stabilizes. It has been found that 
about five to six iterations are generally sufficient to reach the stable 


conditions. 


The program has a provision for specifying the minimum and maximum 


thickness of each element by specifying the proper values of BMIN and 
BMAX in the input Card Sets 11 and 12. The minimum number of laminae in 
a layer is defined by THMIN in the input Card Set 7. The elements can 
be linked to have the same sizes by defining the number of sets and the 
element numbers in each set in the input Card Sets 13 and 14. 

When the structure is to be analyzed only, set LSTCCL = 0 in the 
input Card Set 3. The elements which violate the specified criteria 
can be separated by assigning LCHEK = 1 in the input Card Set 4. This 
control will print out all the elements which violate the strength 
criteria and also the ratio of the actual stress to the maximum allowable 
stress for those elements. If the specified criteria is the maximum 
strain criteria, it will be the es of strains and for the Hill and 
Norris criteria, the effective stress ratio. 

In the output the user has the option to print the stresses in the 
local element coordinate or the fiber coordinate system in a layer (See 
Figure 2), the strains in the fiber coordinate system, the effective 
stress ratios in each layer as defined by the Hill or Norris criteria, 
the average stresses in the 0° fiber coordinate system, and the nodal 
forces in the local or global coordinate system. The control parameters 
for this data are defined in the input Card Set 4. The reactions at the 
boundaries and the sum of the forces and moments about the three coordi- 
nate axes through the origin are printed at the end. This-information is 


generally useful for an equilibrium check. 
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The program listed in the report is limited to solving problems 
with a maximum value of the different variables as given below 


Number of members - MEMBS - 120 


c 


Number of quadrilaterals and triangles - NC - 50 
Number of shear panels and posts - NI - 70 

Number of nodes - JOINTS - 50 

Number of boundaries ~ NB - 50 

Number of loading conditions - LOADS - 2 

Number of layers in the composite elements - NZ - 4 


Number of sets of the materials for shear panels and bars - 
NMAT1 - 4 


Number of different materials used for quadrilaterals and 
triangles — NMAT2 - 4 


Number of elements in the stiffness matrix - NMAX - 5000 
Number of sets of linked elements - LNK - 5 


Number of elements in each linked sets - LMAX - 10 

Core required on the Cyber-7400 - NOS/BE operating system - 114,K 

On the Cyber-7400 a much larger problem can be solved by increasing 
the dimensions of the different variables as indicated by the comments 
in the benginee: 

A short description of the subroutines used in the program is 


given below. 


Program MAIN 

The MAIN program reads the input data, generates the stiffness 
matrix of the whole structure and determines the displacement vector. 
It then uses the information received from the MEMB subroutine to 
calculate the weight of the structure and to modify the design variables 
by using the recurrence relations. At the end it prints out the number 
of laminas in the layered elements, the cross-sectional areas of the 


bar elements and the thicknesses of the remaining elements. 
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SUBROUTINE MEMB 


This subroutine is the longest subroutine and is called twice from 
the MAIN program. It calculates the stresses, strains, and strain energies 
in all the elements. The scaling parameter A used to satisfy the strength 
criteria in all the elements is determined in this subroutine. The membrane 
elements which are susceptible to local buckling are modified here. The 
stresses, strains, nodal forces, etc. are printed out by this subroutine 
after the last iteration. 

The following subroutines are called by MAIN and MEMB to carry out 
various calculations. 

SUBROUTINE SUGD 

This subroutine calculates the scaling parameter for the membrane 
element to satisfy the maximum stress criteria or the maximum strain 
criteria. 

SUBROUTINE AVECT 

This subroutiné defines the Avraatlan of the fibers in the layer of 
a membrane element by a set of components in the local coordinate system 
of the element. 

SUBROUTINE SURFACE 

This subroutine calculates the scaling parameter for a shear panel 
and a bar element to satisfy the maximum allowable stress. 
SUBROUTINE COORD 

This subroutine establishes the local coordinate system of the 
elements. It also determines the transformation matrices needed to 
transform the stiffness matrix and the displacement vector from one 


coordinate system to another. 
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SUBROUTINE ELAS 

This subroutine determines the elastic matrix of the layer in the 
membrane element in the local coordinate system. It also calculates the 
transformation matrix to get the stresses and strains in the fiber 
coordinate system from the local coordinate system. 
SUBROUTINE STRAIN 

This subroutine evaluates the strains in a membrane element from 
the nodal displacements in the local coordinate system. 
SUBROUTINE COMP 

This subroutine calculates the stiffness matrix of a triangular 
membrane element in the local coordinate system. 
SUBROUTINE ELSTIF 

This subroutine determines the stiffness matrix of a bar element 
in the local and the global coordinate system. 
SUBROUTINE ASEMBL 

This subroutine assembles the stiffness matrix for the whole struc- 
ture in the global coordinate system 
SUBROUTINE CONDNS 

This subroutine eliminates the internal dummy node point of the 
quadrilateral and condenses the stiffness matrix from (10x10) to (8x8). 
SUBROUTINE SUM 

This subroutine assembles the stiffness matrix (10x10) of the 
quadrilateral from the stiffness matrices of the four triangles. 
SUBROUTINE CHANGE 

This subroutine which is called from the CONDNS subroutine rear- 


ranges the stiffness matrix of the quadrilateral and the shear panel. 


SUBROUTINE TRNSF 

This subroutine is used to transform the stiffness matrix from the 
local coordinate system to the global coordinate system. 
SUBROUTINE POP 


This subroutine determines the exact core required to store the 


stiffness matrix of the whole structure. It also defines the location 


of the elements on the main diagonal and the first nonzero element of 
each row. 
SUBROUTINE GAUSS 

This subroutine solves the linear algebraic equations and determines 
the displacement vector. 
SUBROUTINE BOUND 

This subroutine eliminates the rows and columns corresponding to 
the 6 undary conditions from the stiffness matrix of the whole structure. 
SUBROUTINE PRNTDR 

This subroutine is called to print out the node numbers, coordinates, 
forces, and the displacements. 
SUBROUTINE REDUCE 

This subroutine eliminates the rows from the force matrix corres- 
ponding to the boundary conditions. 
SUBROUTINE RESTOR 

This subroutine restores the boundary conditions in the force and 
the displacement matrices. 
SUBROUTINE ELSTRS 

This subroutine determines the stress and the strain energy in a 


bar element. 


SUBROUTINE SCOMP 

This subroutine computes the stiffness matrix of a shear panel. 
SUBROUTINE SSRS 

This subroutine determines the stress and the strain energy in a 
shear panel. 
SUBROUTINE BFORCE 

This subroutine calculates the average stress on the the membrane 
element. 
SUBROUTINE SBUCKL 

This subroutine determines the increase in thickness of the mem- 
brane element needed to prevent local buckling of the plate. 
SUBROUTINE ANORM 

This subroutine normalizes the relative design vector. 
SUBROUTINE MULT 

This subroutine is used to multiply the stiffness matrix and the 
displacement matrix to determine the applied forces and the reactions 
at the supports. 
SUBROUTINE ENGS 

This subroutine determines the strain energy in a layer of a 
membrane element. 
SUBROUTINE STRESS 


This subroutine determines the stresses in the local coordinate 


system and the material coordinate system and the strains in a membrane 


element. 
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This subrout ‘ eterming the iling parameter in a membrane 


element for the tre t riteria of Hill or Norri: 
BROUTINE AVG 
Thi subroutine determines the average stresses for a quadri- 


iteral from the stresses in the four triangles. 


he f all the external forces and the moments about the three 
axe massing through the origin are calculated in this subroutine tx 


heck the overail equilibrium. 


This subroutine calculates the average strain energy density for 
all the elements which are specified to have the same thickness for 


linking of variables. 
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SECTLON IV 


PROGRAM INPUT ENSTRUCTIONS 


This section describes the input data required for optimization of 


a structure. 


Card Set 1 PROBLEMS (15) 


Col 1-5 NSTAR Number of problems 


Note: Subsequent set of cards are repeated depending on the number 
of problems to be solved. 


Card Set 2 HEADING (8A10) 


Col 1-80 i) ee An alphanumeric description of the 
problem to be solved. 


Card Set 3 GENERAL PROBLEM RELATED DATA (1615) 
Col 1-5 MEMBS Number of elements 
6-10 NC Number of membrane elements 


(quadrilaterals and triangles) 


11-15 NI Number of shear panels and 
posts (bar elements) 


16-20 JOINTS Number of nodes 


21-25 NB Number of restrained degrees 
of freedom 


26-30 LOADS Number of loading conditions 


31-35 MM Number of degrees of freedom 
per node 


[2 - two-dimensional problem 


3 - three-dimensional problem 


36-40 NZ Number of layers in the 
membrane elements. Each 
layer corresponds to a 
fiber orientation. 


- < . ae : ea eae ae. 


41-45 


46-50 


51-55 


56-60 


61-65 


S 66-70 


INCHES 


KIPS 


NMAT1 


NMAT2 


NFIBER 


NCRTIA 


INCHES = 


KIPS = 


NFIBER = 


NCRTIA = 


30 


Control for units of 
distance 


1 - for inches 


0 - for feet 


Control for the units of 
the applied forces 


1 - for kips 


0 - for pounds 


Number of different sets 

of materials used for shear 
panels and bars (See 
Section IIT) 


Number of different materials 
used for membrane elements 
(See Section IIT) 


Control for the definition 
of 0° fibers 


0 - Projection specified 


11 - Angle 9 is specified 
for each element 


fone 


Control for the strength 
criteria 


1 - Maximum stress (See 
Eq. 6) 


2 - Maximum strain (See 
Eq. 7) 


3 - Hill(See Eq. 8) (NEF in 
Card Set 4 must be 1) 

4 - Norris(See Eq. 9) (NEF in 
Card Set 4 must be 2?) 


Pa 


71-75 


76-80 


1-5 


6-10 


11-15 


NENG 


LINK 


MCONST 


NSTART 


LSTCCL 


LINK 


MCONST 


NSTART = 


31 


Control for the recurrence 
relation 


- Strain energy density 
(See Eq. 24) 


~ Ratio of strain energy 
to maximum strain 
energy density (See 
Eq. 27) 


- Effective stress ratio 
(See Eq. 29) NCRTIA must 
be 3 


- Effective stress ratio 
(See Eq. 31) NCRTIA must 
a be 4 


Control for linkage of 
elements (See Section III) 


0 - No linkage 


1 - linkage 


Control for specified sizes 
for a set of elements 


O - Small number of specified 
sizes (less than 20% of 
the total number of 
elements) 


1 - Large number of 
specified sizes 


Control to read the element 
sizes for the initial design 


0 - Not specified 


1 - Specified 


Number of iterations for the 
stress constraint design. If 
LSTCCL = 0 the given structure 
is to be analyzed only. 


16-20 LLSTCL Number of iterations to 
modify the design variables 
by using the recurrence 
relation (See Section IIT) 
LLESTCL < LSTCCL 


21-25 LNSB Nuinber of iterations for 
modification of the membrane 
elements to avoid local 
buckling (See Section IIT). 


26-30 NBLNCE Control to equalize the 
thickness of + 6 layers 
(last two layers). 


0 - Not required 


NBLNCE = 
1 - Required 


Card Set 4 OUTPUT CONTROL (1615) 


Note: This set controls the nature of the output desired. The program 
will write only that information which the user asks. Most of 
these controls are for the stresses in the quadrilaterals and 


triangles. 
Col 
1-5 NAREA The relative sizes of the 
elements for intermediate 
iterations 
0 - Not required 
NAREA = 
1 - Required 
| 
6-10 NLO The stresses in the layers 
| of the membrane elements in 
the local element coordinates. 
0 - Not required 
NLO = 
1 - Required 
| 
11-15 NZEO The stresses in the layers 
of the membrane elements in 
the material coordinates of 
the 0° layer 
0 - Not required 
NZEO = 
1 = Required 


| 
' 32 


@ 16-20 


21-25 


26-30 


31-35 


36-40 


NSTRN 


NAVG 


NFI 


NEF 


NFOR 


- 


NAVG 


NFI 


NEF 


NFOR 
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The strains in the layers 
of the membrane elements in 
fiber direction 


0 - Not required 


es - Required 
The average stresses for 
the laminate in the 0° 


fiber direction 


0 - Not required 


|i - Required 

The stresses in the layers 

of the membrane elements in 
the individual layer material 
coordinates 


0 - Not required 


1 - Required 


The effective stresses in the 
layers of membrane elements 


0 - Not Required 
1 - Required (See Eq. 8) 
[2 - Required (See Eq. 9) 


The nodal forces in the 
elements 


0 = Not required 


1 - Required (local element 
coordinates) 


'2 - Required (global 
coordinates) 


-_ ee 


41~45 LCHEK Identification of the 
elements which do not 
satisfy the strength 
criteria, when the structure 
is to be analyzed only 
(LSTCCL = 0). 


0 - Not required 


LCHEK = 
a Required 


Card Set 5 DATA TO PREVENT MEMBRANE ELEMENTS 
FROM LOCAL BUCKLING (8E10.4) 


Note: Card set 5 is needed only if LNSB defined in card set 3 is not 
equal to zero. 


PECT (1) Fractions to control the 
fF Doa.ca NS increase in thickness of 
each layer to prevent local 
buckling (See Section IIT) 


NZ 
a PECE (CE) = 1. 


Card Set 6 MATERIAL PROPERTY DATA FOR SHEAR 
sane PANELS AND BARS (4F20.4) 


Note: Units of the allowable stresses in this set should be consistent 
with that of the loads. There will be two cards for each value 
of I. I is the material property number. 


Col 


1-20 E(1) Elastic modulus in psi/10° 


for the Ith material 


21-40 PMU (1) Poisson's ratio for the 
Ith material 


41-60 GSHE (1) Shear modulus in psi/10° 
for shear panels for the 
Ith material. If this input 
is zero, then 


GSHE(T) = a 


2* (1+PMU (1) ) 
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61-80 SSTRT (1) Allowable stress in tension 
| for bars for the Ith material 


1-20 SSTRC(T) Allowable stress in compression 
for bars for the Ith material 


21-40 SSTRS (1) Allowable stress in the shear 
panels for the Ith material 


41-60 SPWTT (1) Density in lbs/in> for the 
Ith material 


I = 1,...,NMATI 


Card Set 7 MATERIAL PROPERTY DATA FOR 
MEMBRANE ELEMENTS (8E10.4) 


Note: The total number of data cards in this set is equal to NZ * NMAT2. 
There will be NMAT2 subset of cards. The Jth subset contains NZ 
cards giving the properties of each layer in sequence. J defines 
the material property number for the laminate. 


Col 


1-10 ee Re Elastic modulus parallel to 
the fiber direction for the 
Ith layer of the Jth material 
in psi/10°, 


11-20 E22(1,J) Elastic modulus transverse to 
the fiber direction for the 
Ith layer of the Jth material 


in psi/10°. 


21-30 ANUL(T,J) Poisson's ratio V0 for the 
Ith layer of the “Jth 


material. 


31-40 ANU2 (I,J) Poisson's ratio Yq for the 
Ith layer of the Jth 
material. If this input is 


zero, then 


£22(1,J) 
E11(1,J) 


ANU2(I,J) = ANUL(I,J) 


41-50 GSH(T,J) Shear modulus for the Ith 
layer ef the Jth material in 
| psi/10.. If this input is 
zero, then 


GSH(1,J) = ——~1LG,J) 
P (2(1+ANU1(1,J)) ) 


} 
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51-60 SPWTC(I,J) Density in lbs/in™ for the 
Ith layer of the Jth 
material. 


61-70 THEK(1,J) Thickness of a lamina of 
the Ith layer of the Jth 
material. This value is 
used to calculate the number 
of laminas in a layer for 
the final design. 


71-80 THMIN(I,J) Minimum thickness of the 
Ith layer of the Jth 
material. This value is 
equal to the minimum number 
of laminas multiplied by 
THEK(1I,J). 


ron kee 
J = 1,...,NMAT2 


Card Set 8 ALLOWABLE STRENGTHS FOR THE 
MEMBRANE ELEMENTS (5Fi5.2) 


Note: The number of cards in this set is equal to NZ * NMAT2 Thus there 
will be NMAT2 subset of cards. The Jth subset contains NZ cards 
giving the allowables of each layer in sequence. J is the material 
property number for the laminate. Units of the allowable stresses 
in this set should be consistent with that of the loads. If 
NCRTIA = 2 (i.e., Maximum strain criteria) this set will contain 
allowable strains instead of allowable stresses. 


Col 


1-15 SSMAX(1,K,J) Tension allowable parallel 
to the fiber direction in 
the Kth layer of the Jth 
material. 


16-30 SSMAX(2,K,J) Compression allowable parallel 
to the fiber direction in the 
Kth layer of the Jth material. 


31-45 SSMAX (3,K,J) Tension allowable transverse 


to the fiber direction in the 
Kth layer of the Jth material. 
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46-60 SSMAX(4,K,J) Compression allowable 
transverse to the fiber 
direction in the Kth layer 
of the Jth material. 


61-75 SSMAX(5,K, J) Shear allowable transverse 
to the fiber direction in 
the Kth layer of the Jth 
material. 


J = 1,...,NMAT2 


Card Set 9 DATA FOR DEFINITION OF THE 0° 
FIBER (3E10.4) 


Note: Coordinates of three points are given. This defines the plane 
of projection and the projection of the 0° fibers on the plane. 
(See Section II(5)). If MM=2, only XD(I) and YD(I) are input. 
If NIFBER=1 in card set 1, this card set is not read. 


Col 
1-10 XD(T) X coordinate of the Ith 
point. 
11-20 YD(T) Y coordinate of the Ith 
point. 
21-30 ZD (1) Z coordinate of the Ith 
point. 
Eh in ea 
Card Set 10 FIBER DIRECTION ANGLES IN LAYERS 
(8E10.4) 
ANGLE(T) Angle in degrees for the NZ 
layers. The first angle is 
Ors 


ae 


Card Set ll DATA FOR MEMBRANE ELEMENTS 


(QUADRILATERALS OR TRIANGLES) 
(15, 13, 12, 415, 5E10.4) 


Note: There is one card for each membrane element (See Figure 4). The 
elements are numbered in sequence but can be arranged in any order 
within this card set. All membrane elements are first numbered, 
then shear panels and bars are numbered. 


Col 
1-5 I Element number 
6-8 KTYPE (1) Type of element 
3 ~ triangle 
KTYPE(I) = 


4 ~- quadrilateral 


9-10 NMAT (1) Material property number 
of the element 

11-15 MA(T) First node number of the 
element 

16-20 MB(1) Second node number of the 
element 

21-25 MC(T) Third node number of the 
element 

26-30 MD(T) Fourth node number of the 
element. For a triangle 
MD(I) = 0 

31-40 BMAX (TI) Maximum thickness allowed 


for the element. If there 
is no limit on size input 
a large number say 100. 


41-50 BMIN (1) Minimum thickness allowed 
for the element. This 
should be consistent with 
THMIN(I,J) in card set 7. 


Note: AWIDE and BWIDE should be consistent with the element supports 
for local buckling. 
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51-60 AWIDE(1) Effective width of the Ith 
element transverse to the 
0° fibers for local buckling 
Ny (compression) critical. 


61-70 BWIDE(1) Effective width of the Ith 
element parallel to the 0° 
fibers for local buckling. 
Ny (compression) critical. 


Note: If NFIBER 


0 in card set 3, then ZANGLE(I) is zero. 


71-80 ZANGLE (1) Angle in degrees which the 
0° fibers make with a line 
joining MA(I) and MB(I). 
See Fig. l(a). 


Te Bo. 


Card Set 12 DATA FOR SHEAR PANELS AND BARS 
(15, 13, 12, 415, 2E10.4) 


Note: There is one card for each element. The element numbering is 
sequential but they can be arranged in any order within this card 
set. These elements are numbered after all the membrane elements 
are numbered. 


Col 
1-5 I Element number 
6-8 KTYPE(1) Type of element 
5 ~ shear panel 


KTYPE(I) = 
2 - bar 


9-10 NMAT (1) Material property number 
of the element 


11-15 MA(T1) First node number of the 
element 
16-20 MB(1) Second node number of the 
element 
( 21-25 MC (1) Third node number of the 
) | element. 


MC(I) = O for bar. 


26-30 MD (1) Fourth node number of the 
element. 
MD(I) = O for bar. 


| 31-40 BMAX (1) Maximum thickness allowed 

for the element. For a 

| bar element BMAX(1) is the 
cross-sectional area. If 

there is no limit on size, 
input a large number, say 

100. 


41-50 BMIN (1) Minimum thickness allowed 
for the element. For a bar 
element BMIN(I) is the 
cross-sectional area. 


b= (NG 2) <.45MEMBS 


Note: Card sets 13, 14, and 15 are read only if LINK = 1 in card set 3. 


! Card Set pS) DATA FOR LINKAGE (315) 
' 
Col 
1-5 LNK Total number of sets of 
linked elements. 
6-10 NSKIN Number of sets of linked elements 
amongst the membrane elements. 
11-15 NINT Number of sets of linked 
elements amongst the shear 
panels and bars. 
Card Set 14 NUMBER OF ELEMENTS LINKED PER SET 
(1615) 
NLINK (TI) Number of elements in each 


linked set. 
I = 1,.«+ sLNK 


Card Set 55. LINKED ELEMENTS (1615) 
NELEM(K, J) Element numbers in each 
linked set. Each subset 
K = 1,...,LNK js started on a new card. 
J = 1,...,NLINK(K) 


Note: Card sets 16 and 17 are read only if LSTCCL = 0 or NSTART = 1 
in card set 3. 
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Card Set 16 THICKNESSES OF THE QUADRILATERALS 
AND TRIANGLES (8E10.4) 
AL(1;L) Thickness of membrane elements. 
Te Le: e ogee Each subset contains one thick- 
Ly Demcwgate ness of the layer for all 


elements. The number of subset 
is equal to the number of layer: 
Each subset is started on a new 
card. 


Card Set 17 THICKNESSES OF SHEAR PANELS AND 
CROSS-SECTIONAL AREAS OF BARS (8E10.4) 


A(1) Thickness of shear panels and 
I = (NC + 1),...,MEMBS cross-sectional areas of bars. 
Card Set 18 COORDINATES OF NODES (15, 3E10.4) 


Note: The node numbering must be sequential but may be arranged in any 
order within this card set. 


Col 
1-5 I Node number 
6-15 X(1) X coordinate of the Ith node 
16-25 Y(T) Y coordinate of the Ith node 
26-35 ZL) Z coordinate of the Ith node 
L ® Lyexe gd OLNio 
Card Set 19 BOUNDARY CONDITIONS (1615) 
IBND (1) Degree of freedom numbers of 
RW Lee vang lee those nodes which are restrainec 
For node k the degree of freedon 
numbers are 3*K-2, 3*K-l, and 
3*K (X, Y and Z respectively) 
for MM = 3 and 2*K-1 and 2*K 
(X and Y respectively) for 
MM = 2. 
Card Set 20 NUMBER OF FORCES IN EACH LOADING 
CONDITION (1615) 
NJLODS (1) Number of load components 
LE @ Lye «5 LOADS in the Ith loading condition. 
41 
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Card Set 21 LOCATION AND APPLIED LOADS 
3(E10.4, 215) 


Note: The number of subsets is equal to the number of loading conditions 


(LOADS). Each subset is started on a new card. 
Col 
1-10 TFR(J) Value of the Jth load. 
11-15 IM(J) Direction of the load 
1 - X direction 
IMCS): = 2 - Y direction 


3 - Z direction 


16-20 JM (J) Number of the node where the 


load is applied. 


J Lycee gNILODS(2) 
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SECTION V 


PROGRAM OUTPUT DETAILS 


This section describes the output from the program. ‘The items are 
discussed in the order they appear. 

1) The input data is printed out in the same format as the READ 
statements in the program. 

2) The output from the POP subroutine concerning the distribution 
of elements in the stiffness matrix is printed. This information is 

a) GROSS POPULATION - total number of elements in the upper 
triangle of the matrix. 

b) APPARENT POPULATION ~ actual number of elements considered 
as non-zero by a given solution scheme. The apparent population repre- 
sents the number of storage locations required for the stiffness matrix. 
If the apparent population is greater than NMAX, the program will write the 
message INSUFFICIENT CORE TO STORE STIFFNESS MATRIX, DIMENSION OF SK 
SHOULD BE (APPARENT POPULATION NUMBER). 

c) STARTING ROW NUMBER FOR EACH COLUMN - The number of the 
row where the first non-zero element occurs in each column. The variable 
defined in the program is IC(1I), I=l, NN(total degrees of freedom). 

d) NUMBERS OF DIAGONAL ELEMENTS IN SINGLE ARRAY STIFFNESS 
MATRIX - For each column I the actual number of elements, ID(1) in the 


upper triangular matrix up to and including that colum, i.e 


ID(1) = Lat Cd b, 


where Bs is the row number given for Column I in (c). Thus for the last 
column ID(LAST) = APPARENT POPULATION. 
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3) TIME USED IN SECONDS - Time the program has spent in execution 
*, 
up to that point. This item is repeated after each iteration, 
4) RELATIVE AREAS OF MEMBERS - These are printed only if NAREA 
in the Input Card Set 4 is l. 

i) The first set includes the relative total areas of all 
members. For composite elements this will be the relative total thick- 
ness of all layers; for shear panels, the total thickness, and for 
bars, cross-sectional areas. The variable defined in the program is 
AUI), t= ty.<s,MEMBS. 

b) The remaining sets include the relative thickness of each 
layer in the quadrilateral and triangular elements. The variable in 
the program is AL(I,J), I = 1,...,NC 

J = 10.03% 


5) WEIGHT-SKIN 


Weight of all quadrilateral and triangular 
elements 


WEIGHT-S-PANELS - Weight of all shear panels 
WEIGHT-POSTS ~ Weight of all bar elements 
LOTAL~WEIGHT - Total weight of the structure. This is 


the sum of above three weights. 


ELE SKIN = Element amongst the quadrilaterals and 


triangles requiring the largest value of the 


scaling parameter to Satisfy the strength 
criteria. 


ABASE - Scaling parameter for the element 
ELE SKIN. 


ELE STRUCT - Element amongst shear panels and bars 
requiring largest value of the scaling 
parameter to sa'isfy the maximum 
allowable stress. 


CBASE ~ Sealing parameter for the element 
ELE STRUCT. 
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STRUCTURE NO - Problem number 

NO OF LOADS - Number of loading conditions 

CYCLE NO - Iteration number 

EQ CYCLE NO - Number of iterations used to equalize 


the scaling parameter in the elements. 
(See Section IIT) 


NSTBLTY - Number of iterations used to modify 
the elements critical in local buckling. 


The data in item 5 is repeated after each iteration. 
6) The following information is given for the minimum weight design 


after completing the specified number of iterations. 


L = element number 

BASEA = scaling parameter 

L=-C = critical loading condition for the element 

CRI = the critical stress in that element 
1 - stress in fiber direction 
2 - stress in transverse direction 
3 - shear stress 
In the case of the maximum strain criteria this 
number will correspond to the critical strain. 
For Hill and Norris criteria CRI = l. 

LAR = layer number where the stress is critical 
PRCNT = percentage of laminae in the different fiber 


directions. 
In the case of shear panels and bars only L, BASEA, L-C are printed. 
7) The following information is given if the structure is to be 
analyzed only. The output is optional. It is printed only if LCHEK = 1 
in the input Card Set 4. 


L = element number not satisfying the strength criteria. 
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ratio of actual state of stress to the maximum 
allowable stress. In the case of maximum strain 
criteria, BASEA is the ratio of strains and for 
the Hill and Norris criteria BASEA is the ratio 
of effective stress. 


L-C = same as in Item 6 
CRI = same as in Item 6 
LAR = same as in Item 6 


In the case of shear panels and bars only L, BASEA, L-C are 
printed. 

8) MINIMUM WEIGHT CYCLE ~ the iteration number which gave the 
lowest weight design. 

9) STRESSES - The following information is given under this 


heading for the minimum weight design 


a) MEMB = element number 

b) NODES = node numbers defining the element 

c) AREA - area of the element (in.”) 

d) THICK - total thickness of the element (in.) 


| 


e) LAYER(THICK) individual layer thickness (in.) 


f) The following output are optional. They are printed only 
if proper control numbers are input in Card Set 4. 


1) STRESSES IN INDIVIDUAL LAYERS (LOCAL COORDINATES) 
2) STRESSES IN INDIVIDUAL LAYERS (ZERO FIBER DIRECTION) 


3) STRAINS IN INDIVIDUAL LAYERS (FIBER DIRECTION) 


, 4) AVERAGE STRESSES (ZERO FIBER DIRECTION) 
5) STRESSES IN INDIVIDUAL LAYERS (FIBER DIRECTION) 


6) EFFECTIVE STRESS 


7) NODAL FORCES 


For items 1 through 6 the stresses or strains are given for each 
loading condition on one line. The first number is the loading condition, 
then the layer numbers, and the stresses 01» Sys Fp are printed out for 
each layer. The format can accommodate a maximum of four layers. In 

the case of average stresses only three stresses are given. The nodal 
forces are printed for each loading condition on one line. When NFOR = 2, 
Fx, Fy, Fz are given in the global coordinates and when NFOR = 1, Fx and 
Fy are given in the local coordinates. The sequence of nodal forces is 
the same as the sequence of numbers used to define the element. 

@in the case of shear panels MEMB, NODES, AREA,and THICK are given 
and for bar elements MEMB, NODES, LENGTH, and AREA are printed out. In 
these two cases since there is only one pertinent stress in the member, 
the loading conditions and stresses are printed on one line under the 
heading of stresses. 

10) DISPLACEMENTS - The following information is given for each ode 
under this heading for the minimum weight design 

a) JOINT - Node number 
|=) ee a re - xX, y, 2 coordinates of the node 
c) FORCE-X, FORCE-Y, FORCE-Z - Applied forces in the x, y, z 
directions. 
d) DISPL-X, DISPL~Y, DISPL-Z - Displacement in the x, y, and z 
directions. 
The quantities (c) and (d) are given for each loading condition. 

11) NUMBER OF LAMINAE IN COMPOSITE ELEMENTS - The number of laminae 

in all composite elements are given. The element number is given in the 


parenthesis before the laminae numbers. 


: LET Pe 


a 12) TOTAL THICKNESS OF ELEMENTS - The actual thicknesses of all 
the elements corresponding to the minimum weight design are listed. 

13) If LNSB > O items 3 to 12 are repeated for the iterations used 
to modify the quadrilaterals and triangles to prevent local buckling. 

14) REACTIONS - The format for this item is the same as Item 10, 
except that under the heading FORCE-X, FORCE-Y, FORCE-Z, the reactions 
are given corresponding to the fixed boundary conditions. 

15) SUMMATION OF FORCES 

a) Applied loads 
LOAD - loading condition Co 


MX - moment of applied forces about x axis through 
the origin 


MY - moment of applied forces about y axis through 
the origin 


MZ - moment of applied forces about z axis through 
the origin 


b) REACTIONS - Same as (a) above except that the quantities 


are calculated for the reactions at the boundaries. 
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SECTION VI 


TLLLUSTRATIVE PROBLEM 


A cantilever box-beam shown in Figure 5 is designed to illustrate 
the use of the computer program. The top and bottom skins consist of 
four layers with fibers in 0°, 90°, +45 and -45°. The coordinates of 


the node numbers shown in Figure 5 are given in the computer output in 


the Appendix. There are 0 nodes and the nodes 37 through )O are 
ce fixed. The total number of elements are 54. There are 18 quadrilaterals, 
18 shear panels ,and 18 posts. The 6 quadrilaterals at the tip of the 
box-beam consist of four layers of graphite epoxy. The 6 quadrilaterals 
in the middle consist of one layer of boron epoxy and three layers of 
graphite epoxy, and the remaining quadrilaterals consist of four layers 
of boron epoxy. The shear panels and posts have different elastic 
properties. 
Number of elements - MEMBS = 54 
Number of quadrilaterals - NC = 18 
Number of shear panels and bars - NI = 36 
Number of nodes - JOINTS = 40 
Number of restrained degrees of freedom - NB = 12 
Number of loading conditions - LOADS = 2 
Number of degrees of freedom per node —- MM = 3 
Number of layers — NZ = 4 
Coordinates of nodes given in feet —- INCHES = 0 
Applied loads given in kips - KIPS = 1 
Number of materials used in shear panels and bars - NMATI1 = 2 
Number of sets of materials used in quadrilaterals - NMAT2 = 3 
Projection of 0° fibers is specified - NFIBER = 0 
Maximum stress criteria - NCRTIA = 1 
j Recurrence relation used (Eq. 24) - NENG = 1 
4 No linkage of elements - LINK = 0 
) No specified sizes of elements - MCONST = 0 
Initial sizes of the e’ements not specified - NSTART = 0 
Total number of iterations - LTSCCL - 15 
Number of iterations to modify the design variables by using 
the recurrence relation - LLSTCL - 10 
Number of iterations to modity the element sizes to prevent 
local buckling - LNSB - 8 
Equal thickness of layers in +45° and -45° required - NBLNCE - 1 
; 49 


if 


Relative sizes of elements required - NAREA = 1 

Stresses in fiber direction needed - NFI = 1 

Nodal forces in global coordinates required - NFOR = 2 

Only the thickness of +45° and -45° layers are to be increased to 
prevent local buckling. Hence PECT(1) = 0.0; PECT(2) = 0.0; PECT(3) = 
O.5; PECECa) = 0.5: 

In Card Set 12, material property number NMATI1 for all bars is 1] and 


for all shear panels is 2. The elastic properties for the two materials 


are as follows 


Material 1 Material 2 
6 ; 6 . 
EB 10.5 x 10° psi 30.0 x 10 psi 
0.3 05.3 
G - an 
2 ; 2 
Allowable Stress in 25.0 kips/in™ 10.0 kips/in 
Tension 
2 2 
Allowable Stress in 20.0 kips/in 10.0 kips/in 
Compression 
. ‘ : : Ae pia 
Allowable Stress in 15.0 kips/in 8.0 kips/in 
Shear 
3 2 
Density 0.1 1lbs/in 0.28 lbs/in 


In Card Set 11 material property number NMAT2 for quadrilaterals 
13 through 18 is 1; quadrilaterals 7 through 12 is 2 and quadrilaterals 


1 through 6 is 3 


Material 1 Material 2 Material 1 
0° layer graphite-epoxy boron-epoxy boron-epoxy 
90° layer graphite-epoxy graphite-epoxy boron-epoxy 
+45° layer graphite-epoxy graphite-epoxy boron-epoxy 
-45° layer graphite-epoxy graphite-epoxy boron-epoxy 
50 
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Elastic properties and allowable stresses for graphite epoxy and boron 


epoxy are as follows? 


graphite-epoxy boron-epoxy 
Buy 18.5 x 10° psi 2a) x 10° psi 
Es, 1.6 x 10° pat 3.5 % 10° ped 
49 0.208 0.25 
Yo4 0.0203 - 
- 6 , 6 - 
G 0.65 x 10° psi 0.93 x 10° psi 
density ASS, Ibs/in> 0725 lbs/in”? 
thickness of laminae -0052 -0052 
Minimum thickness ~0104 -0104 
of a layer 
Allowable stress 
fiber direction “i Fe 
(tension) 139.0 kips/in~ 166.0 kips/in* 
fiber direction 3 7 
(compression) 92.4 kips/in* 86.0 kips/in™ 
Allowable stress transverse 
fiber direction 2 % 
(tension) 4.95 kips/in 6.0 kips/in™ 
fiber direction 5 s 
(compression) 29.7 kips/in~ 11.86 kips/in® 
: 9 
Allowable stress shear 4.68 kips/in” 3.95 kips/in™ 


The projection of 0° fibers is defined by two points (0.0, 0.0, 0.0), 
(1.0, 0.0, 0.0). The third point defining the plane is (0.0, 1.0, 0.0). 
The line joining the first and second points is perpendicular to the line 
joining the first and the third points. 
The four fiber orientations are 0°, 90°, +49$ and -45°. 
The element connections and coordinates of the nodes are given in 
H the output of the program. The sequence of defining the elements is as 


Z shown in Figure 4. 
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The minimum thickness of all elements is .01l inch. 


it 


For all quadrilaterals AWIDE 
For all quadrilaterals BWIDE 


9.0 in 
12.0 in 


tl 


Nodes 37 through 40 are fixed. Hence, the restrained degrees of 


freedom are 109 through 120. 


The first loading condition consists of a 1 kip load at nodes 1, 2 


, 


3 and 4 acting in the negative z direction. The second loading condition 


consists of a 0.5 kip load at nodes 2 and 4 acting in the negative z 


direction. 
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READ DATA 


ASSEMBLE STIFFNESS MATRIX 
OF THE STRUCTURE 


CALCULATE DISPLACEMENTS 


CALCULATE STRAINS, STRESSES 
AND SCALING PARAMETER 
TO SATISFY STRENGTH CRITERIA 


CALCULATE WEIGHT OF 
STRUCTURE 


MODIFY RELATIVE DESIGN VECTOR 
1) RECURRENCE RELATION 

OR 
2) LOCAL STABILITY 


IF NUMBER OF ITERATIONS 
COMPLETED, STOP, OTHERWISE 
GO TO BLOCK 2 


FIGURE 3. FLOW DIAGRAM 
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